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Abstract. We report on a 20 ksec XMM observation of the distant cluster RX J1120. 1+4318 , discovered at z = 0.6 
in the SHARC survey. The cluster has a regular spherical morphology, suggesting it is in a relaxed state. The 
combined fit of the EPIC/MOS&pn camera gives a cluster mean temperature of kT = 5.3 ± 0.5 keV with an iron 
abundance of 0.47 ± 0.19. The temperature profile, measured for the first time at such a redshift, is consistent 
with an isothermal atmosphere up to half the virial radius. The surface brightness profile, measured nearly up 
to the virial radius, is well fitted by a /3-model , with (3 = 0.781q'q4 an d a core radius of 9 C = 0.44+QQ4 arcmin. 
We compared the properties of RX J1120. 1+4318 with the properties of nearby clusters for two cosmological 
models: an Einstein - de Sitter Universe and a flat low density Universe with Qo — 0.3. For both models, the 
scaled emission measure profile beyond the core, the gas mass fraction and luminosity are consistent with the 
expectations of the self-similar model of cluster formation, although a slightly better agreement is obtained for 
a low density Universe. There is no evidence of a central cooling flow, in spite of the apparent relaxed state of 
the cluster. This is consistent with its estimated cooling time, larger than the age of the Universe at the cluster 
redshift. The entropy profile shows a flat core with a central entropy of ~ 140 keV cm 2 , remarkably similar to 
the entropy floor observed in nearby clusters, and a rising profile beyond typically 0.1 virial radius. Implications 
of our results, in terms of non-gravitational physics in cluster formation, are discussed. 

Key words. Galaxies: clusters - Intergalactic medium - Cosmology: observations - Cosmology: dark matter - 
Cosmological parameters - X-rays: galaxies: clusters 



1. Introduction 

Galaxy clusters occupy a unique position in the scenario of 
hierarchical structure formation, as they are still forming 
today. The statistical properties of this evolving popula- 
tion (e.g. the mass distribution function and the corre- 
lations between physical quantities, at various redshifts) 
therefore provide unique constraints on cosmological sce- 
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In the simplest model of structure formation, purely 
based on gravitation, galaxy clusters constitute a homol- 
ogous family. Clusters are self-similar in shape, and pre- 
dictable scaling laws relate each physical property to the 
cluster total mass (or t emperature ) and redshift z (e.g. 
Bryan & Norman |l998t Eke et al. 



1998|) 



From observa- 
tions with ROSAT and ASCA satellites, it is now well 
established that this simple model fails to explain all the 
observed structural and scaling properties of the nearby 
cluster population (Tozzi & Norman 2001 and reference 
therein). The evolution with redshift of these properties 
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is an essential piece of information, still largely missing, 
to reconstruct the physics of the formation process. This 
information is also important for accurate f2o estimates 
based on the evolution of the cluster mass function (e.g. 
Oukbir & Blanchard 1997). To compare the results of flux- 
limited X-ray surveys with the prediction of the various 
theoretical models, it is necessary to understand the re- 
lation between observed quantities, like the cluster tem- 
perature and luminosity, and the virial mass. Precise Slo 
estimate also requires a good understanding of the survey 
selection function, which further dep ends s ensitively on 
cluster morphology (e.g. Adami et al. 200C ). 



Following the pioneering work of the Einstein medium 
Sensitivity Survey (Gioia et al. 1990), several large, 



well controlled, X-ray samples of distant clusters have 
been assembled in the last years, using ROSAT obser- 
vations (see Gioia 200C for a review). The exceptional 
sensitivity of XMM, associated with good spectroscopic 
and imaging capabilities, now allows the detailed analy- 
sis of these clusters, down to what was only a detection 
limit with ROSAT. We present here the XMM observa- 
tion of RX J1120.1+4318, a cluster at z = 0.6 detected 
in the bright Serendipitous High-rcdshif t Arc hival ROSAT 
Cluster (SHARC) survey (Romer et al. |2000|) . This obser- 



vation was made in the framework of the XMM-Newton 
Vt project, a systematic XMM Guaranteed Time follow-up 
of the most distant (z > 0.45) SHARC clusters (Bartlett 
et al. |200l| ). 

The paper is organized as follows. In Sec. ^, we describe 
the data analysis performed to derive the surface bright- 
ness profile, global temperature and temperature profiles, 
which are presented in Sec. |[ In Sec. [|, we compare the 
physical properties of RX J1120. 1+4318 with the predic- 
tions of the self-similar model of cluster formation. The 
scaled emission measure and temperature profiles, the gas 
mass fraction and Lx~T relation are compared to those of 
nearby clusters. This study is made for two cosmological 
models: an Einstein-de Sitter Universe (EdS, £!o = 1) and 
a flat low density Universe (ACDM ) with = 0.3 and 
A = 0.7. In Sec. ||, we study the thermodynamical state 
of this cluster - possible presence of cooling gas, entropy 
content - to further assess the role of non-gravitational 
processes in cluster formation. Section ^| contains our con- 
clusions. 

Unless otherwise stated all errors on the cluster param- 
eters are at the 90% confidence level. A Hubble constant 
of Ho = 50 km/s/Mpc is assumed. At the cluster redshift, 
1 arcmin corresponds to 456 kpc and 561 kpc for the EdS 
model and the ACDM model, respectively. 

2. Data Analysis 

2.1. Data preparation 

RX J1120. 1+4318 was observed for ~ 20 ksec on May 
8, 2001 with the EPIC/MOS and pn camera (using the 
THIN optic al blocking filt er) in Full Frame mode (Turner 
et al. 2001 ; Struder et al. 2001 ) . We generated calibrated 



event files using the tasks emchain and epchain of the SAS 
V5.1. 

We discarded the data corresponding to the periods of 
high backg round induced by solar flares (e.g. see Arnaud 
et al. 2001 ). We extracted the light curves in the energy 



band [10 - 12] keV and [12 - 14] keV for the MOS and pn 
data respectively. In these energy bands, the effective area 
of XMM is negligible and the emission is dominated by the 
particle induced background. We removed all frames cor- 
responding to a count rate greater than 15 ct/100 s (MOS 
data) and 22 ct/100 s (pn data). After this selection, the 
remaining exposure time is 17.6 ksec, 17.9 ksec and 14.2 
ksec for the MOS1, MOS2 and pn observations. 

2.2. Vignetting correction 

The effective area of the XMM mirrors is a function of off- 
axis angle and this vignetting effect depends on energy. An 
additional vignetting effect is due to the RGA obscuration 
for the MOS camera. The vignetting calibration data were 
those available at time of release of SAS V5.1. 

To correct for vignetti ng eff ects, we used the method 
proposed by Arnaud et al. ( 2001 ) for spectra, which can be 
generalized in a straightforward way to profiles or images 
(e.g. Majerowicz et al. 2002 ). For each event, we computed 
the corresponding weight coefficient, defined as the ratio 
of the effective area at the photon position and energy to 
the central effective area at that energy. When extracting 
spectra, image or surface brightness profile, each event is 
weighted by this coefficient]]. These 'corrected' products 
correspond to those we would obtain if the detector re- 
sponse were flat. The on-axis response can then be used 
for spectral fitting and to estimate the physical param- 
eters of the cluster. We used the spectral response files 
ml_thinlv9q20t5r6_all_15.rsp, m2_thinlv9q20t5r6_all_15.rsp 
for EPOC/MOSl,2 and epn_ff20_sY9_thin.rmf V6.1 for 
EPIC/pn respectively^. Those files include both the ef- 
fective area and the redistribution matrix. 



2.3. Background subtraction 

After cleaning for flare events, the XMM background is 
dominated by the Cosmic X-ray background (CXB) and 
the non X-ray background (NXB) induced by high en- 
ergy particles. The former component, which dominates 
at low energy, depends on the observed position in the 
sky and is vignetted by the optics. The latter component, 
which dominates at high energy, is not vignetted but varies 
slightly across the detector and with time. 

To subtract the total XMM background, we used EPIC 
blank sky event files (one for each camera) obtained by 
combining several high galactic latitude pointings. The 
data are cleaned for backgro und fl ares and bright point 
sources are excluded (Lumb [2002 ). The sky coordinates 



1 This vignetting correction can now be done with task 
EVIGWEIGHT of SAS V5.2 



see jittp: / /xmm.vilspa.esa.es/calibration/ 
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Fig. 1. Combined EPIC/MOSl&MOS2&pn image of RX 
J1120 in the [0.3 — 5] keV energy band (linear intensity). 

in the event files were modified using the aspect solution 
of the RX J1120. 1+4318 observation, so that extraction 
can be done in sky coordinates, while insuring that the 
same detector region is considered for both the blank field 
and the source observation. The background level, esti- 
mated from the total count rate in the whole FOV in the 
high energy bands defined above, was found to be 0.85 
times smaller for the source observation than for the blank 
sky data. This is ty pical of the variations observed since 
launch (Lumb 2002), which are explainable by variations 
in Cosmic Ray shielding as the magnetosphere is pumped 
up and down by solar activity. The background subtrac- 
tion for each source product (spectrum or profile) is done 
in two steps. The method is fully described in Appendix [A]. 
We first subtract the corresponding blank field product, 
obtained using the same spatial and energy selection, and 
normalized by the 0.85 factor defined above. For consis- 
tency, the blank field products were obtained using the 
same vignetting correction method as that used for the 
source. The correction factor is thus formally wrong for 
the NXB component. However, as we subtract blank field 
and source count rates estimated in the same region in de- 
tector coordinates this correction factor is the same and 
does not introduce bias. This first step thus allows us to 
subtract properly the NXB. However, we are left with a 
residual CXB component, which is the difference between 
the CXB in the source region and in the blank fields (mul- 
tiplied by the normalization factor above). This residual is 
corrected for vignetting and is expected to be uniform all 
over the field of view. In a second step, we thus subtract 
this residual component, using data in the outer part of 
the FOV, outside the cluster region. 



o 
'—i 



S 10" 

g 

(50 

■ r- 1 

o 

S-H 

3 

on 




0.1 1 
Radius (arcmin) 

Fig. 2. Combined MOS1, MOS2 and pn surface brightness 
profile RX J1120 in the [0.3 - 3] keV energy band. Dotted 
line: Best fit /3-model . Solid line: same convolved with 
the XMM/PSF and binned as the data. Dashed line: On 
axis XMM/PSF, normalized to the central intensity. 

3. Results 

The vignetting corrected image in the [0.3 — 5]keV en- 
ergy band is displayed in Fig. [l|. The data of all EPIC 
cameras are combined. The cluster has a regular spherical 
morphology, suggesting it is in a relaxed state. 

3.1. Surface brightness profile 

We extracted the surface brightness profile of the clus- 
ter in the [0.3 — 3] keV energy band. This band was cho- 
sen to optimize the signal-to-noise (S/N) ratio. We cut 
out serendipitous sources in the field of view and binned 
the photons into concentric annuli with a width of 3.3" 
(3 pixels of the MOS camera) centered on the maximum 
of the X-ray emission for each camera. Since the cluster 
is regular, this center corresponds to the centroid of the 
emission. The three profiles were then summed. The vi- 
gnetting correction and background subtraction was per- 
formed as described above. After subtraction of the cor- 
responding blank field profile, only a contribution from 
the CXB remains (see Appendix |A|). Since the data are 
corrected for vignetting effects, the CXB surface bright- 
ness profile should be constant with radius. The profile 
was indeed found to be flat beyond 4', where we can thus 
consider that the cluster emission is negligible. The resid- 
ual background was thus estimated from the data within 
4' — 7' and subtracted from the profile. Starting from the 
central annulus, we re-binned the data in adjacent annulii 
so that i) at least a S/N ratio of 3cr is reached after back- 
ground subtraction and ii) the width of the bin increases 
with radius, with A(0) > 0.19. Such a logarithmic radial 



4 



Arnaud et al.: XMM-Newton observation of the distant(z = 0.6) galaxy cluster RX J1120. 1+4318 



binning insures a S/N ratio in each bin roughly constant 
in the outer part of the profile, when the background can 
still be neglected. 

The resulting surface brightness profile, S(6), is shown 
in Fig. |[ The cluster emission is significantly detected up 
to i?dct = 3' or 1.37 Mpc for a critical density Universe. 
Beyond that radius it was not possible to create S/N > 3 
annulus of any width. The total count rate within the R < 
3' region and in the considered energy band is 0.30± 0.006 
ct/s. 

We fitted S(8) with a /3-model convolved with the 
XMM PSF, and binned as the observed profile. The PSF 
of each camera is modeled by a normalized King profile, 
with parameters depending on energy and off-axis angle 
(Ghizzardi |2001| , Griffiths & Saxton ^002[ ). The overall 
PSF at each radius is obtained by summing the PSF of 
each camera, estimated at an energy of 1 keV, weighted by 
the respective cluster count rate in the [0.3 — 3] keV energy 
band. The on-axis overall PSF is plotted in Fig. ^ (dashed 
line). The convolution with the PSF takes into account 
the (small) PSF variation across the region considered. 
However, in practice only the MOS PSF variations are 
taken into account. We used the on-axis pn PSF, due to 
the lack of available parametrical fit of the off-axis data. 

The fit of the cluster profile gives /3 = 0.78io;o4 and 



a core radius of C — 0.44 



/+0.06' 
-0.04' ' 



The reduced x 2 is ~ 1 



(x 2 — 20 for 21 d.o.f). The convolved best fit model is plot- 
ted in Fig. U (solid line), together with the corresponding 
(unbinncd and non-convolved) /3-model (dotted line). As 
can be seen, the PSF mostly affects the core of the profile. 

We then estimated the cluster surface brightness pro- 
file corrected for the effect of the PSF, that we will use 
in the following (Sec. |]). In principle, we should have de- 
convolved the observed profile. However since the cluster 
profile is well fitted by a /3-model , we used a simpler pro- 
cedure. For each radial bin, we estimated the ratio of the 
surface brightness corresponding to the non-convolved and 
convolved best fit /3-model . The observed surface bright- 
ness profile was then corrected for the effect of the PSF, 
by simply multiplying the observed S(8) value in each bin 
by the model ratio obtained for that bin. 

3.2. Mean temperature 

The overall MOS1, MOS2 and pn spectra, extracted from 
the event file, are shown in Fig. ^. The spectra are 
corrected for vignetting and background subtracted (see 



Sect. 2.3). To optimize the S/N ratio, the integration re- 
gion for the cluster was restricted to 2' from the clus- 
ter center, and the residual CXB spectrum was estimated 
from the 4' < 6 < 11' region. The spectra are binned so 
that the S/N ratio is greater than 3 a in each energy bin 
after background subtraction. 

The spectra are jointly fitted with XSPEC us ing a red- 
shifte d MEKAL m odel f Mewe et al. |l985|Jl986| ; Kaastra 
Licdahl et al. 1995 ). We let the relative normaliza- 
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Fig. 3. XMM spectra of the cluster from the 9 < 2' region. 
Black (red) [green] points: EPIC/MOSl(2)[pn] data. The 
EPIC spectra are background subtracted and corrected 
for vignetting as described in Sec. |[ Solid lines: best fit 
isothermal model with kT = 5.3 keV, an abundance of 
0.47 times the solar value. 



a common temperature and abundance. When letting the 
hydrogen column density Nn and redshift z free, we ob- 
tained Ah = 2.2 ± 1.2 x 10 20 cm' 2 , in agreement with 
the 21 cm value (A^h = 2.26 x 10 20 cnT 2 from Dickey 
1990|) and z = 0.61 ± 0.03, consistent with 



& Lockman 

the optical value (z — 0.60). We then fixed these param- 
eters to the 21 cm and optical values. The best fit gives 
kT = 5.3 ± 0.5 keV and an abundance of 0.47 ± 0.19. The 
reduced X 2 is ~ 0.94 ( X 2 = 258 for 276 d.o.f). 

The temperature values estimated separately from the 
MOS and pn spectra are consistent within the error bars 
(kT = 5.8± ;? keV for MOS data and kT = 4.5±J^ keV 
for pn data). We note, however that a lower temperature 
is obtained with the pn data than with the MOS data. 
This is likely to be partly due to the remaining calibration 
systematic uncertainties in the XMM spectral responses. 
Indeed, it is known that the EPIC-pn and MOS cameras 
show a relative flux difference of 4% at low energies which 
increases with energy above 4.5 keV, resulting in a MOS 
spectral slope flatter than the pn (Saxton [2002 ; Griffiths 
et al. [2002D. 



tion between the various instruments be free but assumed 



3.3. Temperature profile 

We then extracted the spectra in three concentric annuli, 
centered on the cluster X-ray emission peak and fitted the 
data as described above. The corresponding temperature 
profile is shown in Fig. |J. It appears flat up to 2', within 
the error bars. 

The emission from the central annulus is affected by 
the PSF and may contaminate the other annulii, blurring 
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Fig. 4. Radial temperature profile as a function of angu- 
lar radius, derived from XMM/EPIC data. The horizontal 
line corresponds to the mean value derived from fitting the 
overall spectrum of the region within 2' in radius. 



out gradients if they exist. Moreover the energy depen- 
dence of the PSF, if not taken into account, might in prin- 
ciple bias the temperature estimate, since the photon re- 
distribution is energy dependent. However, for the MOS1 
instrument for instance, the Encircled Energy Fraction 
(EEF) within 15" , the size of the first bin, is already about 
70% (Ghizzardi [200% The flux in the first bin is about 
1/3 of the second bin flux. We thus expect a contamina- 
tion of only ~ 10% in the second bin and we probably do 
not significantly underestimate possible gradients in the 
central part. Furthermore, the EEF varies only by 4.7% 
between 0.3keV and 5 keV, the minimum and maximum 
energy for the spectra. Neglecting the PSF energy depen- 
dence results in a negligible bias (~ 0.5%) at high ener- 
gies. We emphasize that the XMM telescope PSF has a 
very weak energy dependence, so that indeed the correc- 
tion is tiny and we are not facing the problem experienced 
in spectro-imagery with ASCA, for example. 



4. Scaling Properties 

4.1. The self-similar model 

The self-similar model is based on simple assumptions 
for cluster formation, derived from the top-hat spherical 
collapse model. The virialised part of a cluster, present 
at a given redshift, corresponds to a fixed density con- 
trast as compared to the critical density of the Universe 
at that redshift and the internal shape of clusters of dif- 
ferent masses and z are similar. Self similarity applies to 
both the dark matter component and the hot intra-cluster 
medium, whose mass fraction is assumed to be constant. 
Consequently simple scaling laws relate each physical X— 
ray property, Q, to the cluster total mass (or temperature 
T) and redshift, in the form Q oc A(z)T a . 



This simple model thus makes a definitive prediction 
in terms of the evolution of cluster properties. First, the 
normalization (but not the slope) of the Q-T scaling re- 
lations should evolve with z. In particular, for a critical 
density Universe, the virial radius at fixed T should de- 
crease with zas (z + 1)~ 3 / 2 , while the luminosity and cen- 
tral emission measure increase as (z + l) 3 / 2 and (z + l) 9 / 2 
respectively. Furthermore once expressed in scaled coor- 
dinates^, the radial profile of any physical quantity (e.g 
density, temperature) should be the same at all redshifts. 

In this section we will compare the structural and scal- 
ing properties of RX J1120. 1+4318 to the properties of 
nearby clusters. We will consider the scaled emission mea- 
sure and temperature profiles, the gas mass fraction and 
the Lx~T relation. Results are given for both the EdS 
model and the ACDM model^. The physical properties of 
RX Jl 120. 1+4318 are summarized In Table j| 

4.2. The scaled emission measure profile 



Recently, Arnaud et al. ( 2002| , hereafter A AN) studied 
the surface brightness profiles of a sample of 25 distant 
(z = 0.3 - 0.83) hot (kT > 3.5keV) clusters observed 
with ROSAT, with published temperature from ASCA. 
They found that the scaled profiles of distant clusters 
are perfectly consistent with the average scaled profile 
of nearby clusters, for a flat low density Universe. For 
RX J1120. 1+4318 , we will make the same comparison, 
and we only briefly summarize the method. 



4.2.1. The cluster profile 

We used the surface brig htne ss profile, corrected for the 
effect of the PSF (see Sect. 3T). It is converted to emission 
measure (EM) profile: 



EM(r) oc S(9) (1 + z) 4 /e(T, z) 



(1) 



where e(T, z) is the emissivity in the energy band con- 
sidered, taking into account the instrument response and 
interstellar absorption, and r — 9dj±, where dx is the an- 
gular distance. 

The EM profile is then scaled according to the self- 
similar model. We use the standard scaling relations of 
cluster properties with redshift and temperature, with 
the empirical slope of the M gas -T relation derived by 



Neumann & Arnaud (2001 ). This empirical relation is con- 
sistent with the observed steepening of the Ly^-T relation 
and reduces significantly the scatter in the scaled profile of 
nearby clusters. The physical radius is thus scaled to the 
virial radius, R v , with R v oc A z 1/2 (l + z) 3/2 T 1 / 2 and 

3 The considered quantity is normalized according to the 
scaling relation estimated at the cluster temperature and red- 
shift and the radius is expressed in units of the virial radius 

4 The comparison of distant and nearby cluster properties 
does not depend on the assumed Ho value: scaled and physical 
quantities depend on Ho via a multiplying factor, the same at 
all redshifts. 
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Fig. 5. Comparison between the scaled emission measure profile of RX J1120. 1+4318 (data points) and the mean 
scaled profile of nearby hot clusters (full line, see text) for two cosmological models. The error bars on the scaled 
variables take into account the error on the temperature and on the surface brightness profile. The dotted lines 
correspond to the reference nearby profile, plus or minus the corresponding standard deviation. Left panel: Results for 
a EdS Universe (O = 1, A = 0). Right panel: Results for a ACDM model (fl = 0.3, A = 0.7). 



the emission measures by A^ 2 (1 + z) 9 / 2 T 138 (see AAN 
for details). This scaling depends on the density param- 
eter f2o and the cosmological constant A, via the factor 
A z = (A c (O z , A)ft )/(187r 2 ft z ), where A c (fl z ,A) is the 
cluster density contrast at redshift z and f2 z the corre- 
sponding Universe density parameter. Analytical expres- 
sion of A c (fl z , A) can be found in Bryan & Norman (1998). 
The scaled profiles also depend on these cosmological pa- 
rameters via the angular distance used to convert angular 
radius to physical radius. 

The scaled emission measure profiles of 
RX J1120. 1+4318 derived for the EdS model and 
the ACDM model are displayed in the left and right panel 
of Fig. U, respectively. The corresponding virial radii are 
given m Table @. As in AAN, the normalization of the 
Ry-T relation is taken from the simulation Evrard et al. 



(|1996|) 



4.2.2. Comparison with local data 

For each cosmological model, the scaled profile of 
RX J1120. 1+4318 is compared in Fig. | with the corre- 
sponding average scaled profile of nearby clusters. This 
profile (full line), and the typical dispersion around it 
(dotted lines) is derived, as in AAN, from a sample of 
15 hot nearby clusters observed with ROSAT (Neumann 
& Arnaud [1999[ ) . The comparison is performed at radii be- 
yond ~ 0.1i?v, where nearby clusters were found to obey 
self-similarity^]. 

To be meaningful, this comparison requires that the 
relative error in the calibration of XMM/EPIC and 
ROSAT/PSPC can be neglected. Our previous study of 



A1795 (Arnaud et al. 2001), where we performed a com- 
bined fit of the EPIC/MOS and ROSAT/PSPC spectra, 
indicates that the fluxes derived by the two instruments 
match at the ~ ±5% level. Similar relative calibration 



uncertainties were derived by Snowden (2005) 



5 A large dispersion was observed in the cluster core proper- 
ties, which are dominated by non-gravitational physics. 



For both cosmological models, the scaled profile of 
RX J1120. 1+4318 is consistent with the local reference 
profile, taking into account its intrinsic scatter. The 
RX J1120. 1+4318 data thus reinforce the validity of the 
self-similar model. However, some significant discrepancy 
in shape can be noted. At large radii the cluster profile 
falls off more rapidly than the reference profile. This cor- 
responds to its higher (3 value, (3 — 0.78, as compared 
to the universal value of j3 = 2/3, which fits well the 
reference curve (Neumann & Arnaud 1999). The scaled 
core radius, x c = r c /i?v, on the other hand, is similar to 
the value (x c = 0.12) of the reference profile: we obtain 
x c ~ 0.11 (EdS model) and x c ~ 0.14 (ACDM model). 
This might simply be due to intrinsic scatter in the prop- 
erties of distant clusters. A typical scatter of ~ 20% is 
indeed observed in th e slo pes of nearby cluster profiles 
(Neumann & Arnaud 1999 ). A scatter in cluster proper- 
ties is observed as well in numerical simulations of cluster 
formation (Navarro, Frenk & White 1997 ). However, this 
discrepancy could also indicate a systematic departures 
from the self-similar model considered. With data on a 
single cluster, we obviously cannot distinguish between 
these two possibilities. 

It is also interesting to compare in more detail 
the results obtained for the two cosmological mod- 
els. The relative position of the scaled profile of 
RX J1120. 1+4318 with respect to the reference profile is 
higher (in y axis) for the ACDM model than for the EdS 
model. As discussed by AAN, the derived scaled profiles 
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of distant clusters depend on the cosmological parame- 
ters, mostly via the angular distance d\. The typical de- 
pendence is ex d\ at fixed scaled radius. At z = 0.6, this 
d\ factor is about 85% higher for a ACDM model than 
for an EdS model. The d\ dependence is strictly exact 
for a profile of logarithmic slope —3 (e.g. in the clus- 
ter external region for a /3-model with (3 — 2/3). Due 
to the steepest slope of the RX J1120. 1+4318 profile at 
large radii, the effect of varying the cosmological param- 
eters is slightly smaller, of the order of ~ 60% (Fig. ||). 
The effect is smaller than the intrinsic dispersion in the 
nearby profiles at radius below typically 0.2i?v, and be- 
comes only marginally larger than the ~ ±25% dispersion 
at higher radii. As emphasised by A AN, no definitive con- 
clusion on the cosmological parameters can then be drawn 
from the observation of a single distant cluster. We sim- 
ply note that the profile of RX J1120.1+4318 is in better 
agreement with the local reference profile at large radii 
for a ACDM model than for an EdS model. In the latter 
model, the profile appears to deviate more and more from 
the reference profile above ~ 0.4i?y The better agreement 
observed for the ACDM model can simply be an artifact, 
due to the intrinsic dispersion in the cluster properties, 
or again due to systematic departure from the self-similar 
model. However, it is consistent with the statistical anal- 
ysis of the scaled profiles of distant clusters, performed by 
A AN, which clearly favors a low density Universe. 

Finally, we note the statistical quality of the scaled 
emission measure profile, as compared to ROSAT/ASCA 
data. The profile of RX J1120.1+4318 is traced nearly up 
to the virial radius, further out than the mean profile of 
nearby clusters observed with ROSAT and the errors bars 
on most of the data points are smaller than the typical 
dispersion in the reference profile. Note also the very good 
sampling of the profile. Using the same criteria (at least 
3cr detection in each bin), most of the ROSAT profiles of 
distant clusters (AAN, Figure 4) are much more sparsely 
sampled. 



4.3. The scaled temperature profile 

Cluster temperature profiles can now be measured with 
high precision with XMM. First observations of nearby 
clusters indicate that the profiles are isothermal (within 
typically 10%) up to half the virial radius (Arnaud et al. 
2001 ; Arnaud |200T| ). 

With XMM we can now also study the evolution of the 
temperature profiles. The cluster temperature profile, nor- 
malized to the mean temperature, is plotted in Fig. ^ as 
a function of scaled radius (ACDM model). The temper- 
ature profile of RX J1120 is mapped up to about 0.5 Ry. 
This is the first time that a temperature profile can be 
measured at such high redshift. Although the errors bars 
are large, this distant cluster appears isothermal, as ob- 
served for nearby clusters. This is consistent with the self- 
similarity in shape expected in the simple model of cluster 
formation. 



1.5 
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Fig. 6. Scaled temperature profile. The temperature is 
normalized to the mean value and the radius is expressed 
in units of the virial radius (ACDM model, see Table.|]). 

4.4. The gas mass fraction 
4.4.1. VT method 

We first estimated the virial radius (Ry) and mass (My) 
from the best fit temperature and the theoretical Ry-T 
and My-T scaling relations at the cluster redshift, as de- 



scribed in Sec. 4.1. We recall that these relations corre- 



spond to a fixed density contrast at redshift z and are 
derived using the virial theorem (VT): 



G \im v My 
2~fiV 



/3 T kT 



(2) 



where the normalization factor is assumed to be p\ = 
1.05, from the simulation of Evrard et al. ( 1996 ). In that 
case: 



My = 2.98 10 15 A~ 1/2 (1 + z) 



-3/2 



kT 



3/2 



10 keV J 



M (3) 



The best fit temperature corresponds to the emission 
weighted temperature within 2', which is about half of the 
virial radius (see Table |l|). Since ~ 90% of the emission 
lies within that radius, this temperature is a good esti- 
mate of the overall emission weighted value. Within that 
radius, there is no significant evidence of strong gradients. 
Nevertheless, if the temperature profile decreases beyond 
2' ~ 0.5 Ry, the measured temperature would be an over- 
estimate of the mass weighted temperature (considered in 
the simulations). However, if the shape of the temperature 
profile does not evolve with z, the corresponding bias in 
the VT mass estimate would be the same at all redshifts 
and will not affect the following comparison with nearby 
cluster properties. 

The central hydrogen density, ?ih,o[] and the corre- 
sponding gas mass within the virial radius, M sas , are de- 
rived from the best fit /3-model parameters (central sur- 

6 fiH.o = 0.85 n e ,o, where n Ci o is the central electronic density 
for an ionised plasma with the observed abundance 
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face brightness, [3 and core radius), assuming a cluster ex- 
tent equal to the virial radius. The errors on these quan- 
tities, due to that on the surface brightness profile, are 
derived as in Elbaz et al. (1995). The error due to the un- 
certainty on the temperature, which appears in the cmis- 
sivity factor, is negligible. We emphasize that the estimate 
of the gas mass does not require severe extrapolation of 
the data and is therefore robust. The cluster emission is 
detected virtually up to the virial radius for an EdS model 
(-Rdct ~ Rv) and up to about 0.75 Ry for a ACDM model 
(Table 0). The gas mass within that radius is already 70% 
of the gas mass within Ry. 

The central hydrogen density, gas mass, virial radius, 
virial mass and gas mass fraction, / gas = M gas /My, are 
given in Table for the EdS and ACDM models. The 
uncertainty on this last quantity is dominated by the un- 
certainty on the temperature, through the estimate of the 
virial mass. 

4.4.2. BM method 

We also estimated the total mass using the hydrostatic 
equilibrium (HE) equation and the isothermal /3-model 
(BM method): 



G unipMy 
2i?v 



3/3 {Ry/r c f 
2 l + (i? v /r c ) 2 



kT 



(4) 



where (3 and r c are the slope and core radius parameters 
of the gas distribution (derived Sec. 3.1). 

The corresponding virial mass and radius, correspond- 
ing to the same density contrast as in the VT approach, 
are given in Table 0. Comparing Eq. || and Eq. and ne- 
glecting the terms (r c /Ry) 2 , the ratio of the virial masses 
(or radius) derived from the VT and BM method is: 



Ry(BM) ( 3/3 
Rv(VT) ~ ^2^ 



1/2 



My(BM) ( 3 (3 
M V (VT) ~ \2~]h 



3/2 



(5) 



The VT and BM methods give very similar results (see 
Table 0). As compared to the VT estimate, the virial ra- 
dius and mass, estimated using the BM method, are ~ 5% 
and ~ 15% higher, respectively. The central hydrogen den- 
sity is unchanged. The gas mass within the virial radius 
is about ~ 5% higher (slightly larger integration region) 
and the gas mass fraction is decreased by about 10%. 

4.4.3. Comparison with the local gas mass fraction 

The gas mass fraction of RX J1120. 1+4318 can be com- 
pared to the 9 0% confidence region derived by Arnaud & 
Evrard Ql999[ hereafter AE) for hot (kT > 4 keV) nearby 
clusters using the same total mass estimate methods. AE 
foundQ: / gas = 20.1[-2.5, +3.0]% using the VT method 
and / gas = 21.5[-3.5, +4.5]% with the BM method. For 

7 These values were derived for an EdS model, but the lo- 
cal value is not sensitive to the cosmological parameters. The 
definition of the virial region is also slightly different but close 
enough to have no impact on our conclusion: our definition cor- 



Table 1. Summary of the physical properties of 
RX J1120.1+4318 for an EdS model (fi = 1) and a 
ACDM model (fi = 0.3, A = 0.7). 



Model 




EdS 


ACDM 


r c (kpc) 




203+^ 


250±™ 


Rdot (Mpc) 

n H ,o (10 -3 cm" 3 ) 




1.37 
6 q +0 - 6 

°- y -0.5 


1.68 
fi 2+ - 5 

D - z -0.4 






Virial Theorem 


Ry (Mpc) 




1.40 ± 0.07 


2.28 + 0.11 


My (10 14 Mq) 




5.8 + 0.8 


9.5 + 1.3 


M gas (< R Y ) (10 14 


Mq) 


0.85 + 0.04 


1.85 + 0.12 


/gas (%) 




14.6 ± 2.1 


19.6 + 3.0 






HE + isothermal /3-model 


Ry (Mpc) 




1.47 + 0.08 


2.40 + 0.12 


My (10 14 M ) 




6.1 + 0.9 


11.0 + 1.7 


A/ gas (< Ry) (10 14 


M ) 


0.89 + 0.04 


1.94 + 0.12 


./gas (%) 




13.3 ± 2.1 


17.6 + 3.0 


Lx (10 40 ergs/s) 
C 

Cmod {%) 




1.39 + 0.08 

1 7+0.6 

-OA 

2.0 + 0.3 


2.14 + 0.12 
9 4+0.8 

1.65 + 0.25 



both the VT and BM methods, the value derived for 
RX J1120. 1+4318 is perfectly consistent with the AE lo- 
cal value for a ACDM model, but is significantly lower for 
an EdS model. 

However, the temperature of RX J1120. 1+4318 (~ 
5.3 keV) is lower than the median temperature of the 
AE hot cluster sample (~ 8 keV). We are thus not ex- 
actly comparing clusters of same temperature and there 
is some indication that the gas mass fraction varies with 
T. The empirical slope of the M gas -T relation derived by 



Neumann & Arnaud (2001) corresponds to / gas oc T - 44 , 
if the classical My-T relation is not modified. This vari- 
ation is consistent with the variation, / gas oc y°- 41 ±o i6^ 
found by Mohr et al. (|l999|) using the virial mass esti- 
mate. A variation of / gas with T is also readily appar- 



ent in Fig. 3 of AE. 
a possible bias of *- 



Correcting the AE / gas values for 



(5.3/8) c 



gives a local value of 



/gas = 16.8% (VT approach), now in agreement with the 
RX J1120.1+4318 estimate for both a ACDM model and 
an EdS model, although marginally in the latter case. 
The gas mass fraction has been proposed from some 



time (Pen 1997; Rines et al. 1999) as a novel distance indi- 
cator. The method is based on the assumed constancy of 
this quantity with z, while its estimate from X-ray data 
depen ds on the assumed angular distance. As Mohr et al. 
( |2000| ) pointed out, the region considered to compare / gas 
at different redshifts must be fixed in terms of scaled ra- 
dius, with the evolution of the virial radius properly taken 
into account, to avoid biases due to the variation of the 
gas mass fraction with physical radius. The comparison 
performed above, where the gas mass fractions are esti- 
mated within the virial radius, does not suffer from this 



responds to a density contrast of 178 at z~0 for an EdS model, 
while the value quoted are estimated at a density contrast of 
200 
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bia^]. However, as illustrated above, the comparison is fur- 
ther complicated by possible variation of / gas with cluster 
temperature. A much better understanding of the / ga s _ T 
relation, both in the local and distant Universe is required, 
before any conclusion can be drawn on the cosmological 
parameters. 



4.5. The L x -T relation 
4.5.1. Cluster luminosity 

We computed the bolometric luminosity, Lx, within the 
virial radius (given Table |l|). The observed count rate in 
the [0.3-3] keV band, obtained by integrating the surface 
brightness profile up to the detection radius, is converted 
to bolometric luminosity using the best fit MEKAL spec- 
tral model and the instrument response. The contribution 
beyond the detection region is estimated from the best fit 
/3-model , but is totally negligible (0.7% for the ACDM 
model). The error on Lx includes both the statistical error 
on the count rate and on the temperature. 



4.5.2. Evolution of the Lx~T relation 



As in AAN and in Sadat et al. (199S), we divided this 
luminosity by the luminosity estimated from the local 
Lx~T relation of AE and the cluster temperature: Lx = 
1.15 10 45 (T/6 keV) 2 88 . The resulting factor C is given in 
Table [j] for the two cosmological models considered. The 
contribution to the error due to the uncertainties on T 
and Lx are summed quadratically in the log space. 

This C factor can be compared to the evolution 
of the normalization of the Lx~T relation, C mo d(z) = 
A^ 2 (1 + z) 3 ^ 2 , expected in the self-similar model 
(AAN). This theoretical value, estimated at the cluster 
redshift, is given in Table | for the EdS and ACDM mod- 
els. The 'error' bars corresponds to plus or minus the in- 
trinsic scatter in Lx estimated by AE. The luminosity of 
RX J1120. 1+4318 is in good agreement with the expected 
evolution for an EdS Universe, and in the upper range of 
the expected evolution for a ACDM Universe. For an EdS 
model, the luminosity is also marginally consistent with 
no evolution. 



4.5.3. Comparison with recent works 

Several groups quantified the evolution of the Lx~T re- 



1998 



lation using ROSAT/ASCA dat a (Sad at et al 
Reichart et al. |1999| ; Fair ley et al. |2000[ AAN). There is a 
general consensus that no significant evolution is observed 
for an EdS model. However, as emphasized by AAN, both 
the luminosity estimate and the theoretical evolution de- 
pend on the assumed cosmology and a different conclusion 



8 Note that the dependence of / gaa on flp and A in the present 
approach (given in Neumann & Arnaud (2000)) is indeed dif- 



ferent than the d \ dependence derived by Pen ( 1997 ) and 



is reached for a ACDM model. In that case, AAN found 
a significant evolution, consistent with the self-similar 
model. Our results agree with this finding. They are 
also consistent with the results of Reichart et al. (1999): 

^(0.91-1.12^™ 



2.0±g;| extrapo- 



they found an evolution of C = (1 
in the redshift range z < 0.5 or C 
lated at z = 0.6 for the ACDM model considered here. 
This is similar to the value, C = 2.4+q 8 , derived for 
RX J1120.1+4318 . 

On the other hand, Borgani et al. ( 20011: ) , using recent 
CHANDRA data up to z — 1.26, did not find any evidence 
of significant evolution, for the EdS model but also for the 
ACDM model. We note first the large uncertainties on the 
CHANDRA data, especially above z = 1. Moreover, there 
is a large dispersion, with data points above and below the 
local Lx~T curve for the ACDM model. To illustrate this 
point, let us consider the two clusters RX J0848+4456 
(z = 0.57, kT = 3.6 + 0.5 keV) and MS1137.5+6625 
(z = 0.78, kT = 5.7t° i keV )- The y are both relaxed clus- 
ters, with precise kT measurements, and are in the same 
redshift and temperature range than RX J1120. 1+4318 . 
The C value of MS1137.5+6625, C = 2.0±g;£, is perfectly 
consistent with the expected evolution: C mo d(z) = 1.88 at 
the cluster redshift and is similar to the value obtained for 
RX J1120.1+4318 . On the other hand, the luminosity of 
RX J0848+4456 is particularly low for its measured tem- 
perature: C = 0.6^Q 2, corresponding to a negative evolu- 
tion and of course inconsistent with the expected positive 
evolution of C mo d(z) — 1.6. The measured evolution is 
sensitive to systematic uncertainties on the temperature 
(C* cx T" 3 ). RX J0848+4456 was observed with ACIS-I 
(Holden et al. 



2001), for which there are still large cali- 



Rines et al. (L999) 



bration uncertainties, especially at low energies. This is a 
potential worry for temperature measurements of high z 
clusters with this instrument. The temperature estimate 
of MS 1137.5+6625, observed with ACIS-S (Borgani et al. 
|2001bD , is a priori more secure. MSI 137.5+6625 was also in 
the AAN sample and the Chandra/ACIS-S data are con- 
sistent with the ASCA data. First cross calibration stud- 
ies of ACIS-S and XMM/EPIC also show a go od agree- 
ment between these instruments (Snowden 2002 ) . Further 
cross-calibration studies and accurate measurements on a 
larger cluster sample are definitively required to assess if 
the presently observed dispersion is real. 

4.6. Consistency of the results 

In our comparisons with the nearby cluster properties, we 
found that the scaled EM profile of RX J1120.1+4318 , 
as well as its gas mass fraction and luminosity, are con- 
sistent with the predictions of the self-similar model of 
cluster formation. Since all three observables are not in- 
dependent - they are directly related to the gas density - 
the consistency of our results is a priori not surprising. It 
is, however not entirely trivial, as we discuss now. 

A good match of the scaled EM profile of 
RX J1120. 1+4318 with the local profile, would imply that 
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i) the shape of the gas profile is similar to that of nearby- 
clusters ii) the EM and thus gas density scales with z 
as expected. Therefore its luminosity (integrated emission 
measure profile) and gas mass (integrated density profile) 
would naturally be found to follow the standard evolution. 
Since the standard evolution of the gas mass is the same 
than the evolution assumed to compute the virial mass, we 
would also find a gas mass fraction consistent with the lo- 
cal value. However, the emission measure, luminosity and 
possibly gas mass fraction also depend on the tempera- 
ture. We emphasize that independent studies of the evo- 
lution of these quantities are thus expected to yield consis- 
tent results, but only if one considers a consistent scaling 
with T for the three quantities. The present analysis is 
consistent in that sense. As mentioned above, the theoret- 
ical model used to scale the EM profiles (M gas oc T 194 ) is 
consistent with the slope of the local Lx~T relation (used 
to normalize the cluster luminosity) and corresponds to 
/gas oc T 44 (used to estimate the local gas mass fraction 
at the cluster temperature). 

In spite of the general agreement with the self-similar 
model, some differences appear when comparing the re- 
sults of the EdS and ACDM models. The analysis of the 
gas mass fraction and scaled profiles would rather favor 
a ACDM model. However for such a model, the luminos- 
ity is in the upper range of the expected evolution and 
a better agreement is observed for an EdS model. This 
is actually a direct consequence of the specific shape of 
RX J1120. 1+4318 . The concentration of the gas distri- 
bution of RX Jl 120. 1+4318 is slightly larger than the ref- 
erence profile, yielding a higher luminosity as compared 
to the gas mass. These differences can also be understood 
by looking at the scaled profiles. As outlined above, the 
scaled profile of RX J1120. 1+4318 is in a good agreement 
with the nearby profile at high radii for a ACDM model. 
For an EdS model, the scaled profile lies more and more 
below the reference curve at r <; 0.4 i?v, indicating a lower 
gas density than expected. Since the external regions con- 
tribute most to the gas mass (~ 75% of the mass lies 
beyond 0.3i?v), this explains why the gas mass fraction 
for an EdS model is marginally too low, even as compared 
to the corrected AE local value, while the agreement is 
better for a ACDM model. On the other hand, for an EdS 
model, the cluster scaled profile matches better the local 
profile in the central region (see Fig.^), which contributes 
most to the X-ray luminosity. 

5. Core properties and non-gravitational effects 

5.1. Radiative cooling 

The XMM data suggest that RX J1120.1+4318 does not 
host a cooling flow. There is indeed no indication of a tem- 
perature drop in the central bin and a good fit is obtained 
with a /3-model , down to the center, with no indication 
of central excess. 

We computed the central cooling time, t coo \ cor- 
responding to the measured central density: < coo i = 
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Fig. 7. Entropy profile of RX J1120.1+4318 , derived 
from the best fit /3-model and assuming isothermality. 
The radius is scaled to the virial radius for the EdS (dot- 
ted line) and ACDM (full line) models. The error bars 
corresponds to the error bars on the temperature profile, 
plotted at the central radii of the bins considered. 

2.9 10 10 yrs-v/T/riH, where T is in keV and tih in cm~ 3 
(Sarazin p.986| ). We obtained t cool ~ 1. 10 10 yrs. The age 
of the Universe at the cluster redshift is t = 6.5 10 9 yrs for 
the EdS model and t = 1.1 10 10 yrs for the ACDM model. 
The cooling time, even at the cluster center where it is 
close to its minimum value, is thus similar to or larger 
than the age of the Universe. We thus do not expect a 
strong cooling flow in the center, in spite of the apparent 
relaxed state of the cluster. 



5.2. Entropy and non-gravitational heating 

As mentioned in the introduction, the simplest self-similar 
model of cluster formation fails to explain all the observed 
properties of the nearby cluster population. A definitive 
evidence of breaking of self-similarity is the entropy ex- 
cess (the "entropy floor" ) detected in the core of cool sys- 
tems with a baseline entropy of about 70 — 140 keV cm 2 
(Ponman et al. 



1999; Lloyd-Davies et al. 2000). The ori- 



gin of this break of similarity is not yet understood. 

With the present data, we can for the first time study 
the entropy profile of a z = 0.6 cluster, ie look at cluster 
entropy evolution. Note that RX J1120. 1+4318 appears 
to be in a relaxed state and the entropy distribution is 
not likely to be affected by shocks induced during re- 
cent mergers. The cluster entropy profile, as a function of 
scaled radius, is plotted in Fig. for the EdS and ACDM 
cosmologies. The entropy S = T/n^ 3 is estimated from 
the best fit /3-model and the cluster mean temperature. 
The uncertainty on the profile is dominated by the un- 
certainty on the temperature distribution. Typical errors, 
corresponding to the temperature error in each bin of the 
temperature profile are indicated in the figure. 

One observes a nearly constant entropy level of about 
140 keV cm 2 within the central r < 0.1 i?v region and a 
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rising entropy profile beyond that region. This shape is not 
surprising, the entropy profile shape reflects the shape of 
the inverse of the density profile, which has a core radius of 
about 0.1i?v- The core entropy is, however, in remarkable 
agreement with the entropy "floor" measured in nearby 
clusters. The same agreement was noted by Arabadjis et 
al. ( pOO^ ) in the case of EMSS 1358+6245 at z = 0.33. 

If confirmed on more clusters, this coincidence has po- 
tentially substantial implications for the physics of clus- 
ter formation. We first note that it is consistent with the 
expectation of the early pre-heating scenario, where the 
gas is pre-heated at a given entropy level befo re col lapse 
(e.g. Borgani et al. |2001a| ; Tozzi & Norman |200l|) . At 



the beginning of the cluster formation, the entropy floor 
prevents shock heating and the collapse is adiabatic. The 
initial entropy is thus preserved in the core. When the 
mass of the system increases, the infall of incoming shells 
becomes supersonic and a shock regime begins (gravita- 
tional heating). The external profile follows the classical 
'self-similar' rising profile, as we observe. The core entropy 
of RX Jl 120. 1+4318 suggests that the entropy floor was 
already established at high z: at least z > 0.6 and proba- 
bly much before, since in a hierarchical scenario the core 
of the cluster must have collapsed earlier. Note also that 
the metal abundance found for this distant cluster is sim- 
ilar to the abundance observed in nearby clusters (e.g. de 
Grandi & Molendi 2001 and reference therein). If early 
galactic winds are responsible for the gas pre-heating (e.g. 
Kaiser 1991; Evrard & Henry 1991), we do expect that 



this is accompanied by an early enrichment. Finally, our 
data suggest that pre-heating in hot clusters has mostly 
an effect on the core properties up to high z. This might 
explain why AAN verified self-similarity of the EM profile 
of hot clusters above 0.1i?v up to z=0.8. 

Radiative cooling has been also proposed to explain the 
entropy "floor" : it can remove low entropy gas from near 
the cluster center, tri ggerin g the inflow of higher entropy 
material (Pearceet al. |2000D . RX J1120.1+4318 , contrary 
to most nearby relaxed clusters (and EMSS 1358+6245), 
is a cluster, for which radiative cooling has probably not 
yet affected the core properties (see above) . This suggests 
that radiative cooling cannot be the dominant process in 
the establishment of the entropy floor. 

6. Conclusion 

XMM-Newton data allow us to measure, with un- 
precedented accuracy at a redshift of z = 0.6, the 
gas and temperature distribution of the distant cluster 
RX J1120. 1+4318 . The cluster has a regular spherical 
morphology, suggesting it is in a relaxed state. The cluster 
temperature is kT = 5.3 ± 0.5 keV and there is no signif- 
icant evidence of temperature radial gradient up to half 
the virial radius. The surface brightness profile, measured 
nearly up to the virial radius, is well fitted by a /3-model 
with p = 0.78±g:gS- 

The RX J1120. 1+4318 data reinforce the validity of 
the self-similar cluster formation model. For both an EdS 



and a ACDM model, the scaled emission measure pro- 
file beyond ~ 0.1i?v is consistent with the nearby cluster 
reference profile, taking into account its intrinsic scatter. 
Consistently, the gas mass fraction is in agreement with 
the local value (although marginally for an EdS Universe), 
and the luminosity of the cluster, taking into account its 
temperature, is consistent with the expected evolution of 
the Lx~T relation. 

There is no evidence of a cooling flow at the cluster 
center, in spite of its apparent relaxed state. This is con- 
sistent with the estimated cooling time, larger than the 
age of the Universe at the cluster redshift, indicating that 
radiative cooling has not yet affected the cluster proper- 
ties. The entropy profile shows a flat core with a central 
entropy of ~ 140 keV cm 2 , remarkably similar to the en- 
tropy floor observed in nearby clusters. This favors early 
pre-heating models for the establishment of this entropy 
floor. 

A statistical study of large sample, as the SHARC clus- 
ter sample, covering a wide range of redshift and luminos- 
ity (thus temperature) is required to better constrain the 
physics of cluster formation, from the evolution of clus- 
ter properties. However, the present data already demon- 
strate that XMM provides the statistical quality required 
by such study. 
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Appendix A: Background subtraction method 

In this section we detail the background subtraction method. 
This method applies after screening for flares. We suppose we 
have a 'template' event file, for the background estimate, ob- 
tained by collecting several high latitude observations. In the 
following, the data are supposed to b e corr ected for vignetting 
effect, as described in Arnaud et al. (2001) and in section 2.2: 



each event detected with energy E at location x, y is weighted 
by a coefficient w(x, y, E) which is the ratio of the effective area 
at position (x,y) to the central effective area, for the energy 
E. The Cosmic X-ray Background (CXB) varies across the sky, 
but can be considered as uniform at the scale of 30', the size 
of the field of view. The non X-ray background (NXB) is not 
uniform in the FOV, but is not vignetted. 

For the template file, the corrected count rate, measured 
at a given location (x,y) and energy E, T(x,y,E), is the sum 
of the CXB and NXB contributions: 

T{x, y, E) = TcxbOe, y, E) + T N xb(x, y, E)w{x, y, E) (A.l) 

Since the data are corrected for vignetting effects, the CXB 
component is, apart from statistical fluctuations, uniform over 
the FOV and corresponds to the average CXB for the blank 
field observations: 



T C x B (x,y,E) =T C xb(E) 



(A.2) 
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The NXB component is the NXB for the observations, 
T^xs{x,y,E), multiplied by the weight factor. 

Similarly for the observation data set, the count rate is: 



Orx,y,E) 



with 



S{x,y,E) +OcxB{x,y,E) 
+Onxb(o;, y, E)w(x, y, E) 



C xB(x,y,E) = C xB(E) 



(A.3) 



(A.4) 



where S(x,y,E) is the source contribution and Ocxb(-B) cor- 
responds to the CXB at the pointing position of the consid- 
ered observation. Since the CXB varies across the sky, a priori 
Ocxb(-E) Tcxb(E). The quiescent NXB can be reasonably 
considered to have the same spectral and spatial characteristics 
for all obs ervat ions but there is evidence of long term variations 
(D.Lumb 2002). We can thus write: 



Onxb(i, y,E) = Q Tnxb(», y, E) 



(A.5) 



The normalization factor, Q, is estimated by considering the 
total count rate in the whole FOV in the high energy band, 
where any cosmic X-ray emission is negligible. 

The background subtraction for each source product (spec- 
trum or profile) is done in two steps. In a first step, we subtract 
the corresponding product of the blank field observation, ob- 
tained using the same ener gy a nd s patia l sele ction , no rmal ized 
by the factor Q. From Eq. TO, Eq. O, Eq. O, Eq. Hand 



Eq.A.5 



0(x,y,E)-QT(x,y,E) 



S(x, V ,E) + Ocxb(E) 

-Q Tcxb{E) (A.6) 



and the NXB contribution is thus removed. The remaining 
CXB component, which can be considered uniform, can thus 
be estimated, using data outside the source region, where: 

0(x,y,E)-QT{x,y,E) = Ocxb{E) - Q T C xb{E) (A.7) 

It can then subtracted in a second step. 

To illustrate the method, let us consider that we want to 
estimate the source spectrum in a particular region, Reg: 



S Reg (E) = y S(x,y,E) 



Beg 



(A.8) 



In a first step, we extract the region spectrum from the ob- 
servation data set. The spectrum extracted in the same region 
from the template file is then subtracted, after normalization 
by the factor A n0 rm above. The resulting spectrum is: 

W£) = ^20(x,y,E)-Q^T(x,y,E) (A.9) 

Reg Reg 
= Sr OS (E) 

+A Reg (Ocxb (E) - Q Tcxb (E)) (A. 10) 

where An eg is the surface of the extraction region. We are 
thus left with a residual X-ray background component, which 
is the estimate of the difference between the CXB spectrum 
in the observation and in the template file (multiplied by the 
normalization factor). 



In a second step, we do the same for a region Reg' outside 
the source. The resulting spectrum is: 



Inc B >(E) = A Kcg ,(OcxB(E) - Q Tcxb(E)) 



(A.ll) 



The source spectrum can thus be estimated by subtracting this 
spectrum to the spectrum obtained in Step 1, after norma liza- 
tion to t he size of the extraction region. From Eq. 
Eq. A.ll one get: 



A. 10 



and 



iRcg 



A 



Reg' 



(A.12) 



The same method can be applied to extract the source sur- 
face brightness profile, in a given energy band. We extract the 
surface brightness profile from the observation data set. The 
surface brightness profile extracted in the same energy band 
from the template file is then subtracted, after normalization 
by the factor Q above. The resulting surface brightness profile 
is an estimate of: 



I(r) = 5(r) + (So c 



- Q s Ta 



(A.13) 



where S(r) is the source surface brightness profile and 5o CX b 
and 5t cxb are the CXB surface brightness in the energy band 
considered for the observation and template file respectively. 
The residual CXB contribution can be estimated from the mea- 
sured profile, I(r), in the region outside the source and sub- 
tracted. 
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